Application of Minimal Subtraction Renormalization to Crossover Behavior 
near the '^He Liquid- Vapor Critical Point 
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Parametric expressions are used to calculate the isothermal susceptibility, specific heat, order pa- 
rameter, and correlation length along the critical isochore and coexistence curve from the asymptotic 
region to crossover region. These expressions are based on the minimal-subtraction renormalization 
scheme within the model. Using two adjustable parameters in these expressions, we fit the theory 
globally to recently obtained experimental measurements of isothermal susceptibility and specific 
heat along the critical isochore and coexistence curve, and early measurements of coexistence curve 
and light scattering intensity along the critical isochore of '^He near its liquid-vapor critical point. 
The theory provides good agreement with these experimental measurements within the reduced 
temperature range |i| < 2 x 10~^. 



I. INTRODUCTION 



It is well known that thermodynamic quantities 
exhibit singularities asymptotically close to the crit- 
ical point. The power-law behavior of these singu- 
larities, characterized by critical exponents and the 
concept of universality and scaling, have been suc- 
cessfully described by renormalization-group (RG) 
theory. Away from the asymptotic region, thermo- 
dynamic quantities of real physical systems deviate 
from simple power-law behavior. However, RG the- 
ory can still provide insight in understanding critical 
crossover behavior. 

There are two main field-theoretical 
renormalization-group schemes to treat critical- 
to-classical crossover phenomena. Dohm and 
co-workers developed the minimal-subtraction 
renormalization (MSR) scheme while Bagnuls 
and Bervillier developed the massive renormaliza- 
tion (MR) scheme Q. Both of these theories used 
the Borel resummation technique to describe the 
crossover behavior of the 0"* model in any 0{n) 
universality class and in three dimensions. The 
difference between the two schemes was discussed 
in ref.||l[]. These field-theoretical crossover theories 
were improved over the years as asymptotic theories 
became more accurate ||^. Recently, Larin et al. 
improved the MSR expressions for the specific heat 
and compared their results with the superfluid 
helium [n ~ 2) system [Q. Bagnuls and Bervillier 
have also improved their theory to match the 
recent asymptotic values for exponents and leading 
amplitude ratios Both renormalization schemes 
can provide crossover functional forms for thermal 
properties with a minimal set of fluid-dependent 
adjustable parameters. However, a direct compar- 
ison between these recent theoretical predictions 
and different experimental measurements near the 
liquid- vapor critical point (n = 1) has been lacking. 



In this paper we will present a direct comparison 
between the MSR field theoretical crossover func- 
tions and various experimental measurements near 
the liquid-gas critical point of '^He. The comparison 
using the MR theory will be published elsewhere. 

The paper is divided into two parts. In the first 
part, we briefly summarize the MSR functional ex- 
pressions for susceptibility, specific heat, coexistence 
curve, and correlation length from previous work 
llB H' 0' il- addition, we derived within the 
MSR framework new RG functional expressions for 
the asymptotic critical amplitudes of the suscepti- 
bility and coexistence curve as well as the first coef- 
ficient in a Wegner expansion for susceptibility, spe- 
cific heat, coexistence curve, and correlation length. 
From these expressions, universal amplitude ratios 
for the 0(1), 3-dimensional system are calculated 
and compared with the most recent values from Bag- 
nuls et al. and Fisher et al. 0. 

The second part of the paper includes the re- 
sults of MSR functional fits to experimental mea- 
surements. In our previous work, we analyzed the 
isothermal susceptibility of '^He along the critical iso- 
chore above T,. using theoretical expressions based 
upon the minimal-subtraction scheme [p^ . In this 
work, we combine that analysis with susceptibility 
measurements along the coexistence curve and spe- 
cific heat measurements along the critical isochore 
[pl]| . Measurements of coexistence curve and the 
light scattering intensity near the critical point of 
■^He 1 12, |3| are also analyzed. 



II. THEORETICAL EXPRESSIONS 

The Hamiltonian for the model in three dimen- 
sions [d ~ 3) is 

H^^ j d\{\ro4>l + ^(V0o)' + «o0o} > (1) 
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where 0o is the order parameter field, whose statis- 
tical mean value is the physical order parameter of a 
given system. The parameter uq is the fourth-order 
coupling constant. The parameter rg is related to 
the reduced temperature t = (T — Tc)/Tc by 



To = aot, 



(2) 



where ao is a nonuniversal constant. It is important 
to note that the total Hamiltonian is the sum of 
H — + Hq , where Hq is the analytic background 
free energy. Since the liquid-vapor critical point has 
a single component order parameter and belongs to 
the 0(1) universality class, we have n = 1. 

The dimensionless bare order parameter field (po 
and the bare coupling parameters uq and tq are 
renormalized to [Eqs. (S2.11)and (S2.12) [|] 

^ = Z4u,e)-i/2^o (3) 
u = Zu{u,ey^ Z^{u,e)^ A3U0 (4) 
r = at = Zr{u,e)^^ Got , (5) 

where ^3 — (47r)~^ is a geometric factor and e = 
4 — (i = 1 for dimension d — 3. The Z factors are 
associated with their respective field-theoretic func- 
tions 



Cr{u) = 


ndf,lnZr{u,e) ^ ^ , 


(6) 


CA^) = 


fidf,\nZ^{u,ey^\^ , 


(7) 


I3u{u) = 


u 


'-!+ fid^{Z-^Zl)l 


, (8) 



where the index means differentiation at fixed tq, 
ipQ, and Uq. 

By introducing a flow parameter I, the effective 
coupling u{l) satisfies the flow equation 



du{r) 



(9) 



The flow parameter I is related to the correlation 
length by 



(10) 



with being an arbitrary reference length. The 
flow parameter I = corresponds to the Ising flxed 
point u{l — 0) — u* , which is determined from 
Pu{u*) = 0. The effective coupling r{l) satisfies the 
fiow equation 



' di 



(11) 



The flow parameter / = 1 is an arbitrary reference 
point, at which the nonuniversal initial values are 
u{i = 1) = 7i and r(l) = r ~ at. 



The ficld-thcorctic functions C.riu),^^{u) and 
Pu{u) in Eqs. (||)-(H) are known up to five- loop or- 
der in perturbation expansions around u = 0]. 
However, the expansions do not converge away from 
u = 0. Hence Borel resummations were used on the 
expansions to calculate the values of these functions 
over the range < u < u*. For most investiga- 
tions, see for instance ref. ^, |^, ||], only the function 
values at the fixed point u* were calculated using 
the Borel resummations on the five-loop expansions. 
The function values over the range < u < u* were 
obtained using up to two-loop order expressions with 
extrapolation terms added in order to reproduce the 
values at the fixed point [Q . For a system of dimen- 
sion d = 3 and single component order-parameter 
n = 1, one obtains 



Cr{u) 

C0(w) 



aiu"^ — a2U^ , (12) 
(13) 

-u + mu^{l + aiu)l{l + a5u). (14) 



I2u - 120m 

-24^2 



03 u 



Here oi through 05 are the coefficients for the ex- 
trapolation terms with values listed in Appendix 
Using these functions and the flow equations, ther- 
mal properties along the critical isochore and coexis- 
tence curve can be calculated from the asymptotic to 
crossover regions using the initial values for Eqs. (^) 
and (pi]), u — u{l — 1) and a — ail — 1), and the 
arbitrary length scale (x~^ in Eq. (|lO|). 



A. Reduced temperature 

Within the MSR scheme, the expression for the 
reduced temperature in terms of the fiow parameter 
/ can be derived as follows. The reduced tempera- 
ture t and the fiow parameter / can be linked using 
Eqs. (S4.25), and (S4.26) in ref.§ and Eq_rH2.9) of 
ref . , together with the solution of Eq. (|ll|) 



^(0 = r(l)exp Ic.y 

1 

= a|i|exp / Q-jj- 



2i2 



with 



b+Hl)) 
b-Hl)) 



Qiu{l)), 
l~Q{u{l)). 



(15) 



(16) 
(17) 



Here is for T > and is for T < T^. 
Krause et al. B determined a one-loop expression 
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plus a higher-order extrapolation for Q{u) given by 
[Eq. (K3.5)] 



Q{u) = l + bqu'^ In(cQM) . 



(18) 



where bq and cq are the extrapolation coefhcients 
with the values given in Table [II in Appendix 

By adding and subtracting = Cr(w*) in the inte- 
grand of Eq. and using the identity = 2— C*, 
where v is the critical exponent of correlation length 
5, one arrives at 



/dl' 
(C.-C)7r = ^^±(0M'^'/''- (19) 



By rearrange Eq. ( |T9| ) , one obtains the following ex- 
pression for the reduced temperature 



with 



and 



|t| =6±(/)to/'/"exp[-F,(0], 



to = — exp[^^,(l)], 
a 



Fr{l) 



' 

WC.K)-C.(w*) , , 

au . 



(20) 



(21) 



(22) 



I3uiu') 

B. Susceptibility 

1. General expression 

The following expressions for the dimensionless 
susceptibility Xt = X± were given in refs. |0, |^ 
respectively for T > [Eq. (K2.7)] and T < 
[Eq. (H2.16)], 



X± = 



^iH^f±{u{l)) 



1 ^ 

exp /C0y • (23) 



The amplitude functions, /±, were expressed to 
two-loop order plus a higher-order extrapolation, 
[Eqs. (K3.1) and (H4.2)], to give 

09 

/+[^.] = 1 - -u\l + b^u) (T > T,), (24) 

/_[w] = [l - 18m+ 159.56M2(H-dxu)]~^ (T < T^), 

where b^ and are the extrapolation coefficients 
with the values given in Table [II in Appendix 



The minimal rcnormalization factor, Z^, in 
Eq. (H) is given by [Eq. (K A12)] 



U 

Z(j){u)~^ — exp j du 



Pu{u') 



The expression 

X± = Xo 



j-^h exp[-f0(M(O)] 
f±{uil)) 



(25) 



(26) 



can be obtained by adding and subtracting Q = 
C,^{u*) in the integrand of Eq. (p3|), and using the 
relations Q — — r/ [|l[ and 7 = i^(2 — 77), where 77 is 
the critical exponent of the fluctuation correlation 
at the critical point and 7 is the critical exponent of 
susceptibihty. In Eq. ( p6| ) 

Xo = A*"^ Z^{u) exp[i^0(l)] , 



and 



''^7[C0("(O)-C0K)] 



I 

nil) 



(27) 
(28) 



XoM 



Pu{u') 

i-v), 

2-r, exp[~j^0(M(O)] 



Using Eq. ( |lO[ ) and 7 = iy(2 — 77), Eq. ( P6[) can be 
rewritten as 



(29) 



Thus, in the asymptotic regime (i^0 —>■ 0) 

X±^De~" (30) 
with a nonuniversal proportionality constant D = 

2. Critical Amplitudes 

Within the pure cj)'^ model, the standard Wegner 
expansion for the susceptibility is given by 

x± = r±ir^(i + r± |t|^ + r± \tf^ + ■■■), (31) 

where are the leading asymptotic critical am- 
plitudes, Vf are the first Wegner expansion ampli- 
tudes above and below the transition, and A is the 
correction-to-scaling exponent. The details of the 
derivations of the leading and first Wegner critical 
amplitudes are given in Appendix Here we list 
the derived expressions for the critical amplitudes. 



Xo (bltoV 



f±{u*) 



77— + 



f± 



(32) 



(33) 



with A = h'LU and to — df3u/du\u*. 
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C. Specific Heat 

1. General expressions 

The total specific heat is usually separated as 

C± - + C± , (34) 

where the term Cb > represents an analytic "back- 
ground" contribution from the analytic background 
free energy _ffo, and represents the critical con- 
tribution from order parameter fluctuations. Here 
is for the specific heat above Tc along the criti- 
cal isochore, ' — 'is for below Tc in coexisting phases. 

The critical specific heat derived from the 
Hamiltonian expressed in Eq. (|^) has two representa- 
tions within the MSR scheme. These two represen- 
tations are derived via multiplicative and additive 
renormalization as detailed in ref. ||]. The most re- 
cent work by Larin et aZ.Q used the representation 
via additive renormalization that we will use in this 
paper. 

The critical specific heat per unit volume near 
Tc is expressed by [Eqs. (S2.36) or (L3.3)]|, | 



^0 =^'^''^1^ / ^-^exp-iJ^ 



(35) 



= iaV'^3i^±(^i(0)exp ^ du'. 
The amplitude functions K±{u) are given by 



K±{u) = F±{u) ~ A{u) . 



(36) 



The functions F± (u) for n = 1 can be expressed by 
a two-loop calculation plus a higher-order extrapo- 
lation, [Eqs. (K3.4) § and (H4.4)§] 



F±[u] 



-1 - 6u{l + hpu) (T > Tc) 
{2u)-^ - A{1 + dpu) {T<Tc) 



(37) 



where bp and dp are the extrapolation coefficients 
with the values given in Table |l| in Appendix 0. 
The function A{u) in Eq. (pq) is governed by 



dA{u{l)) 



dl 



AB{u{l)) + {l-2Qr{u{l))}A{u{l)), 

(38) 

with A{u = 0) = -45 (u = 0). The function 
B{u) has been calculated to 0{u^) for any given n 
[Eq. (L2.21)] [Q. However the five-loop Borel resum- 
mation of B{u) was only performed for n = 1 at u*. 
Hence a new extrapolation term with a coefficient 
&B is added to the two-loop expression Q| in order 
to satisfy B{u* , n — 1), 



B{u) 



%l + bBu)u^ 



(39) 



At the fixed point u — u* , ldA{u{l))/dl = 
Pu{u)dA{u)/du = since /3„(u*) = 0, and Eq. (||) 
leads to 



A* = A{u*) 



Ai^B{u*) 



(40) 



where 2(* — 1 = —ajv 0] is used with ol being the 
critical exponent for specific heat at constant vol- 
ume. 

The integral in the exponential of Eq. ( p5| ) can be 
rewritten as 

f (2C. ~l) = 2j'^{Cr- C) + In l-"/" 

(41) 

using 2C* — 1 = —ajv. The expression for the specific 
heat from the additive renormalization can now be 
rewritten as 



cl 



167r/i 



K±{l)r"/'' exp [2Fr{l)] exp [-2T^(1)] 



= Col-"^" exp [2F,,(0] K±{u{l)) , (42) 
where Fr{l) is given by Eq. (^ and Co is defined as 

(43) 



2 3 

Co = exp [-2^,(1)1 - ^ 



167r/i 



2. Critical Amplitudes 

The standard Wegner expansion within the pure 
model for specific heat can be written as 

C± = A^\t\-''il + Af\t\^ + Af\t\^^ + ---) 

+ Bcr + CB (44) 

where Bcr is a constant background induced by long- 
range correlations between the fluctuations. The ex- 
perimentally measured constant background is the 
sum of Bcr and the analytic background, Cb- 

The expression for the Wegner expansion of the 
specific heat via multiplicative renormalization was 
derived and given in Eqs. (S4.23) and (S4.24) of 
ref. §. In Appendix we derive the expressions for 
the critical amplitudes and the critical background, 
Bcr, for the representation via additive renormal- 
ization, using the technique that is consistent with 
the one used for susceptibility. The results of these 
derivations are, 



A^ = Co {bltor {Fl - A* 



(45) 



A* - Fl 



2v 



-{2B' -A*C) 



(0 ^^'r 



w (blto) 



A ' 



(46) 
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- ^ = A{u) -A* + -^{2B' - A*a{u* - u) . 
Go A — a 

^ ^(47) 

The variables with a prime in Eqs. (Hq) and ( p7| ) are 
derivatives with respect to u. The right hand side of 
Eq. ( ^ ) is negative for any given u. Hence one has 
Bcr < since Co > from Eq. (^). 



with the leading critical amplitude and the first Weg- 
ncr amplitude being respectively 

Bo^Mb-tor^if;)'^^ (54) 

^1= Pt P TT^ + T^ 7 Ta- 55) 

V 6- UJ 2/0 2wy (6* to) 



D. Coexistence curve 



E. Correlation length 



1. General expressions 

In the liquid- vapor coexisting phases below Tc, the 
density difference ^Pl,v = Pl,v/ Pc — 1 is the statis- 
tical mean of the order parameter field (</>). There 
is no asymmetry between Api and Apy within the 
(j)^ model. Schloms and Dohm have given an expres- 
sion for the square of the physical order parameter 
[Eq. (S3.10)] i, 

(0)2 = A^Z^{u)U{u{l_))C^ exp C0— ■ (48) 

Here Z^{u) is given in Eq. (p5|). The amplitude func- 
tion f^{u) is expanded in one-loop with an extrapo- 
lation term, [Eq. (H4.1)]0|, to yield 



/^(u) = (8M)"\l + d0M). 



(49) 



The correlation length below Tc is linked to the flow 
parameter by I- = By combining these 

expressions and following the derivation of Eq. (26), 
one has 

(0)2 = 4 &'''U{u{l^)) cM-m-)] , (50) 

where l + r/ = 2/3/z^is used and 

</)2 = (4^)-VZ0(w)exp[F0(l)], (51) 

with (3 being the critical exponent of the order pa- 
rameter. 

2. Critical Amplitudes 

Using Eq. ^ to replace fj^'" in Eq. (|o|) and 
the scaling relations 7 = v{2 — rj), a — 2 — 3i>, and 
a + 2f3 + J = 2, one has 

(0) = ±0o<o''|t|'3[5_(Z_)]-^'[/0(L)]V2 

X exp[/3F,(/_)] exp[-F0(;^)/2] . (52) 

Expanding 6_(Z_), f^{l-), Fr{l-), and F^fL) in the 
same manner as described in Appendix H, one ob- 
tains the Wegner expansion for coexistence curve. 



ApLy = ±So|i|'^(l + Bi \t 



(53) 



Using Eq. ( pO[ ) to express I in terms of the ex- 
pression for dimensionless correlation length (if /i^^ 
is taken dimensionless) is derived from Eq. (nOl) as 

= M"' [b±{u{l))toV exp[-z.F,(0] . (56) 
around u(l) ~ u* to 



An expansion of Eq. (| 
0[{u - u*)2] leads to 



1 + 



(57) 



LU b± 



\t\' 



This equation is identical to Eq. (S4.8) of ref.[p. By 
comparing Eq. ( |57| ) to the standard Wegner expan- 
sion form, 



(58) 



one obtains the leading amplitudes and first Wegner 
correction amplitudes of the correlation length 



^± = i^-'ibitoY , 



C b', 



u 

u* ib' 



±toy 



(59) 
(60) 



F. Universal amplitude ratios 

Even though the leading amplitude and sub- 
sequent Wegner expansion coefficients are fluid- 
dependent, certain combination ratios of these am- 
plitudes are universal. From the equations for the 
first Wegner amplitudes of the specific heat, suscep- 
tibility, coexistence curve, and correlation length, 
one notices that the system-dependent part, {u — 
u*)/{b'^to)'^, is the same in every expression. There- 
fore the ratio of any of these first Wegner amplitudes 
is universal based on the MSR cj)* model. These uni- 
versal ratios have been given for the specific heat in 
ref. [||, [ Eqs. (H), del), and (||)]. In this paper we 
derive the other universal ratios based on the MSR 
0^ model. From Eqs. ( ^ and (^, one has the 
universal amplitude ratios for susceptibility. 



J- n 



/+(«*) 



4.94, 



(61) 



6 



rr 



A 

' I,) 



/+ 



= 0.228. 



^ uj ^ b- ~^ ^ f- 

_ _ (62) 
From Eqs. ( [45|) and (Q), one obtains the universal 
amplitude ratios for specific heat, 



An 



b 



aFl 



A + 

^ = 1.07. 



= 0.535. 



(63) 



(64) 



Use of the scaling relation a + 2/? + 7 = 2 and the 
combination of Eqs. (^2|), (^5|), and (jsj) leads to a 
universal ratio 



(65) 



agoxo(J";-^*)(^;^o)"+^ 
0§(&*_to)-2^/;/; 

a(fe;)"+T(6!L)-2/5(F; - A*) 
4^7* 



= 0.0580. 



From Eqs. ( |59| ) and (|60|), the universal amplitude 
ratios for the correlation length are, 



^= I =1.42. 

Co" V fe- 



et 



A C 


K 




b+ 




b'_ 







= 1.10. 



(66) 



(67) 



From Eqs. (^5|) and (p9D, one has the universal 
relation between the amplitude of specific heat and 
correlation length 



aFl) 



0.0206 
0.0134 



(T > TJ 
(T < T,) 



(68) 



where the scaling relation a + iv — 2 has been used. 
Equation (H) is identical to Eq. (S4.22) of ref.[|. 
The evaluation of the right-hand side uses the con- 
stants given in Appendix A natural extension 
of Eq. (68) is the universal relation between spe- 
cific heat and the correlation length throughout the 
crossover region. Using Eqs. ([ic|), (po|), and (42), 
and the scaling relation a + 3;^ = 2, one has 



167rk|2 



(69) 



TABLE I: The values of various universal amplitude ra- 
tios. The calculation for this work uses the values of the 
am plit ude functions at the fixed point u* given in Ta- 
ble |v| and the values of the critical exponents given by 
Guida and Zirm- Justin j^. 



Amplitude 



ratios 


This work 


B and B" 


F and Z*" 




4.94 


4.79 ±0.10 


4.95 ±0.15 


^0 Mo 


0.535 


0.537 ±0.019 


0.523 ± 0.009 




1.42 




1.89 ±0.015 


Rc 


0.0580 


0.0574 ± 0.0020 


0.0581 ±0.0010 




0.0206 


0.0196 ± 0.0001 


0.0188 ±0.00015 




0.0134 




0.0053 ± 0.00025 




0.228 


0.215 ±0.029 




A+/A- 


1.07 


1.36 ±0.47 






1.10 






Bi/rt 


0.76 


0.40 ± 0.35 





"Bagnuls and Bervillier 
'Fisher and Zinn [Hi 



Since there is no fluid-dependent parameters appear- 
ing on the right hand side of Eq. (|69|) , the product of 
the critical specific heat and the cubic of the corre- 
lation length is universal for any given temperature 
throughout the crossover region. 

Table Q lists the various universal amplitude ra- 
tios derived from the minimal-subtraction renormal- 
ization scheme, Bagnuls and Bervillier's massive- 
renormalization scheme and other methods, such 
as £-expansion, summarized by Fisher and Zinn 
The values given by Bagnuls and Bervillier are 
closely matched to the values given by Guida and 
Zinn- Justin after the readjustment of the Borel 
resummation criteria Q . Noticeable differences exist 
in Table | among various theories. In attempting to 
explain these differences, two factors are important 
to note. First, we are unable to evaluate the uncer- 
tainties of the universal ratios since the uncertainties 
on the Borel resummations at the fixed point u* for 
most of the amplitude functions were not given in 
previous studies. Secondly, Eqs. (m), (M), and (|60| ) 
use the derivatives of Eqs. (|l|), (|l^ (|l§);(|2|), (|3^, 
and ( |39|) which could have sizable systematic uncer- 
tainties. These equations were only obtained from 
two-loop calculations and extrapolated to the five- 
loop fixed point values with adjustable constants. 
Hence it is desirable to have these derivatives cal- 
culated at the fixed point with Borel resummations. 
Then the extrapolation coefficients can be more ac- 
curately reconstructed, leading to the estimates of 
the first Wegner coefficients with less uncertainties. 
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III. FIT TO EXPERIMENTAL 
MEASUREMENTS 

The expressions within the MSR model are para- 
metric for susceptibihty, specific heat, coexistence 
curve, and correlation length along the critical iso- 
chore and coexistence curve. We made a variable 
change of I = exp(— x) in solving those expressions 
numerically, x was discretized with 1000 data points 
over the range of — oo < a; < cxd to obtain the solu- 
tion for u{x) over the range of < u(x) < u* . For 
each thermal property versus reduced temperature, 
say x+ vs t, 1000 data points were calculated for a 
look-up table of x+{xi) vs t{xi) with j = 1, • • -1000. 
The intended property was then obtained for a given 
reduced temperature t using a cubic spline. 

This MSR model has three system dependent pa- 
rameters, u, /i, and a which fix the scales for u{l), 
^(l), and t{l) in Eqs. |, [o|, and ^ Here u{l), ^(0, 
and t{l) are defined implicitly as functions of the RG 
fiow parameter Since I is eliminated in final solu- 
tion, it is clear that one of the three amplitudes is 
redundant. This should not be mistaken as a min- 
imal number of three fitting parameters for a com- 
plete equation-of-state while we only fit the thermal 
properties along the critical isochore and coexistence 
curve. In this paper {/i, a} are chosen as fitting pa- 
rameters for a prefixed u because their combination 
only appears in the amplitude of the parametric ex- 
pressions, such as to, Co, etc. The u value is 
chosen based on the consideration that the expres- 
sions for the first Wegner amplitudes were derived 
by ignoring higher order terms in [u{l) — u*]. There- 
fore an accurate determination of the first Wegner 
amplitudes requires u to be close to u*. 

Besides {/x, a}, the critical temperature can also 
be a fitting parameter. Another adjustable param- 
eter is required for the analytic background contri- 
bution to specific heat. In fitting the experimental 
data, j/expt, to theory, we minimize 



2 / VcxptiXi) 2/thoory (s^i, o) 

i=l 



(70) 



Here a is an array of fitting parameters with the 
standard error a given by 



dy 

dx 



(71) 



The partial derivative in Eq. (71) is evaluated nu- 
merically in each fitting iteration. The x in Eq. ( |70| ) 
is temperature. In our experiment, the sample tem- 
perature was determined from a resistance measure- 
ment having an approximately 10 /^K uncertainty, 
i.e. ax — 1 X 10~^ K. In fitting the measurements of 



isothermal susceptibility and specific heat, we assign 
ay = k X 2//100, assuming a k% uncertainty in the 
measurement. 

The goodness of a fit is characterized by the value 



N -M' 



(72) 



where N is the number of data points and M is the 
number of fitting parameters. 

All the experimentally measured quantities were 
made dimensionless by expressing them in units of 
appropriate combinations of the '^He critical tem- 
perature Tc = 3.315 K, critical density pc = 
0.04145g/cm'^, and critical pressure Pc = 1.14 Pas- 
cal. 

The experimental susceptibility, xt = 
p{dp/dP)T, is scaled by p1/Pc to obtain the 
dimensionless susceptibility, Xt = XtPc/Pc- The 
physical order parameter, Ap = p/ Pc — 1 is already 
dimensionless. The measured heat capacity had 
units of [C] = J/K. It was then divided by the fiuid 
volume to have a unit of [pCy] = J/(cm'^K). Since 
the energy unit is [J] = [P][V], a dimensionless 
specific heat was obtained as Cy = pCyTc/Pc with 
P^/Tc = 0.03463J/(cm3K). 

The critical specific heat per unit volume 
(divided by Boltzmann's constant ks) near Tc is 
given by Eq. (^2|). The volume scaling factor is 
vo = ksTc/Pc, thus the length scaling factor is 



lo = t'o^' = [kBTc/Pc) 



1/3 



(73) 



For ^He, one has lo — 7.36 A. It is assumed that is 
dimensionless in the MSR (f)^ model expressions with 
^0^^ being the scahng factor, one defines ^* = S,/lo- 



A. Fit to susceptibility measurements 

The susceptibility along the critical isochore (p = 
Pc) was determined using PVT measurements from 
both sides around pc- The susceptibility along the 
coexistence curve was also determined using PVT 
measurements. For example, Xt^ obtained from 
PVT measurements for p > p^coex ^'^'^ Xt^ 
tained from p < Pcoex- 

Since Xt varies sharply as p ^ pc and p Pcoex , 
the dominating uncertainty in Xt comes from the 
uncertainty in locating either pc for measurements 
above Tc and p = pcoex for measurements below 
Tc- Above Tc, the inflection point was well con- 
fined by the data from the both sides of pc- Be- 



low Tr 



P 



Pc 



was determined from a kink 



in P versus p curve. However, this kink becomes 
less pronounced as T — > Tc. Based on our ob- 
servation, we assigned the susceptibility uncertain- 
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lO"*" 10'' !()'■* 10-'' 10"- 10"' 10"* 10"- 10"^ 10"' 

\T/T,-l\ l-u/u* 



FIG. 1: Fit of MSR <p'^ model to ^He susceptibility mea- 
surements for both T > Tc and T <Tc. u/u* is fixed to 
be 0.999, -while Tc, ^l, and a were adjusted. The solid line 
is the best fit. The dot-dashed straight lines represent 
the asymptotic predictions from the fit. 



FIG. 2: The susceptibility fitting quality parame- 
ters ^ and a, and the resultant critical amplitudes versus 
fixedu/it*. V'^ and F^ were calculated using Eqs. ( js^ ) 
and (|33[) for each set of {u, fi, a} from the fit. 



ties to be (t^^{T > T^) = 0.02xt{T > T^) and 
<j^^{T <T,) =Q.Ixt{T <T,). 

The results of fitting the susceptibility measure- 
ments for both T > Tc and T < to the MSR 
expression in Eq. ( p6| ) is shown in Fig. |l|. The sus- 
ceptibility was scaled by \t\'^ in order to provide a 
more sensitive representation of the crossover behav- 
ior and the fitting quality. The dot-dashed straight 
lines represent the asymptotic predictions from the 
MSR fit. The uncertainties in the amplitudes were 
deduced from the uncertainties of /i and a in the fit. 

Figure ^ shows and {/^, a} versus (1 — u/u*). 
It is interesting to see that the goodness of the fit 
remains unchanged over a wide range of (1 — u/u*). 
This verifies that only two out of the three fiuid- 
dependent parameters are relevant fitting parame- 
ters. We note here that no improvement in x^ "was 
made within scatters when u/u* was also free to be 
adjusted in the fit. 

Also shown in Fig. ^ are the calculated Tq, and 
Ti versus (1 — u/u*). F^ reaches a plateau for 
(1 — u/u*) < 4 X 10""^ since the analytical expres- 
sion for F^ was derived by ignoring 0((1 — u/u*)'^) 
or higher order terms. Therefore in this work we fix 
the value of u/u* = 0.999. 

Furthermore, /i oc (1 — u/u*)'^'^ and a oc (1 — 
u/u*)^''"^ for (1 - u/u*) < 0.1 implies a strong 
correlation between fx and a when susceptibility is 
chosen in fitting the MSR theory to experimental 
data. The above power-law dependence of /i and a 
on {1 — u/u*) comes from the fact that is nearly 
a constant for (1 — u/u*) < 1 x 10^^. Thus from 
Eq. (H) to ~ ^i^/a cx (1 - u/u*)^/^ w (1 - u/u*)^. 



B. Fit to specific heat measurements 

The specific heat, Cy, near the '^He critical point 
was measured using a heat pulse method. The tem- 
perature change could be measured very accurately 
using a magnetic susceptibility thermometer with 
1 nK resolution. For T > Tc, temperature equili- 
bration was very fast due to the "piston effect" , 
and the uncertainty in the measured Cy "was ~ 1%, 
i.e. ac^{T > T,) = 0.01Cv(r > T,). For T < T„ 
equilibration underwent critical slowing down as the 
fluid approached Tc. The slowing-down was due to 
the mass transfer at the meniscus between liquid 
and vapor. Since the sample cell was not perfectly 
adiabatic due to its mechanical support and electri- 
cal wires, there was some heat loss from the cell to 
the surrounding during the long equilibration. The 
uncertainty in measuring Cy "was typically 5%, i.e. 

In fitting Cv measurements to the MSR (j)* model, 
an additional adjustable parameter, Cb, appears in 
Eq. (|4^). By treating Cb as a constant within a 
small reduced temperature range around Tc, the true 
crossover behavior described by the MSR (j)'^ model 
can be revealed. Figure ^ shows a fit of the Cy 
measurements for both T > Tc and T < Tc to the 
MSR expression in Eq. (p2[). The fit was limited 
to the reduced temperature range \t\ < 2 x 10^^ 
as indicated by an arrow in the figure. The agree- 
ment between the experimental measurements and 
the theory is good. The uncertainties in the criti- 
cal amplitudes and fluctuation-induced background 
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r^. = 3.315668 ± 0.000025 K 
j(/h* = 0.999 

= (1.51 ±0.17) xlO-» 
a = 0.122 ±0.005 

e 2^^i±.....-..Z7. 


A„ = 6.97 ±0.83 
A7=1.1±0.1 
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A|= 1.2±0.1 
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Cg = 3.85 ± 0.05 
B„ = -5.5±0.8 


Xi= 1-193 
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FIG. 3: The dimensionless specific heat at constant vol- 
ume versus reduced temperature. The symbols repre- 
sent the experimental measurements. The solid line is 
the best fit. The dot-dashed straight lines represent the 
asymptotic predictions from the fit. The arrow indicates 
the fitting range |t| < 2 x 10"^. 



were error-propagated from the uncertainties of /i 
and a. 

The fluctuation-induced bac kgrou nd specific heat, 
-Bcr, was calculated from Eq. (C16) in Appendix^. 
Its absolute value is close to that of Cb- As a result, 
the combined background specific heat is close to 
zero for ^He as first demonstrated experimentally 
by Brown and Meyer |15t. 



C. Fit to Coexistence Curve Measurements 

The best data for the ^He coexistence curve was 
compiled in a recent paper by Luijten and Meyer 
[ p2t . We apply the MSR 0* model to the coexistence 
curve using these data as shown in Fig. |^. The fit 
was limited to the range 

The lower bound was so chosen since the measure- 
ments were affected by the gravity effect for |t| < 
6 X 10~^ due to a large cell height (4.3 mm) used 
in that experiment [n6|. The upper bound was so 
chosen since the (/)^ model was developed for critical 
phenomena and did not include analytic behavior 
associated with a system approaching absolute zero 
temperature. The standard deviation for |ApL_y| 
was approximated based on the percentile deviation 
in Fig. 5 of ref. 0, namely 1% at |t| = 6 x 10"^ and 
0.2% at |t| > 1 X 10~^. The standard deviation for 
reduced temperature was (5T = 1 x lO^^K divided 
by Tc = 3.3155K. 

The solid line in the figure represents the best fit 
with only {^,a} adjusted. The predicted Bq — 1-02 
is consistent with the reported Bo — 1-02 in ref. [p2|. 



The agreement between the model calculation and 
the experimental data is satisfactory over the fitting 
range. The systematic difference between the MSR 
(^"^ model calculation and the measurements over the 
fitting range also exists from other theoretical model 
calculations |p^ . 

The systematic deviation between the MSR 0'* 
model calculation and the experimental data over 
the fitting range may be due to the fact that there 
was no proper background contribution included in 
the analysis. We attempt in this paper to include 
the effect of the order parameter saturation as a 
possible background contribution. The saturation 
of order parameter at absolute zero temperature has 
been studied by Povodyrev et al. for an ideal 
Ising model. We propose an empirical expression 
that is consistent with that study for the limiting 
behavior at \t\ ~ 1. Not only does the order param- 
eter saturate to a constant value but its slope also 
approaches zero at \t\ = 1. In the case of the liquid- 
vapor system where the physical order parameter is 
the normalized density difference from the critical 
value, the saturation value is also unity. Our em- 
pirical expression, satisfying this limiting behavior, 
is 

ApL.v = (0)exp(-|t|/ti)±62 [l-exp(-|tl/t2)] • 

_ (74) 
In fitting the expression in Eq. ( |74| ) to the exper- 
imental data, only ti is adjusted while t2 and 62 
are solved for a given ti through the constraints 
^Pl.v ~ ±1 and dAp^y /dt = at \t\ — 1. Near the 
critical point, the exponential damping of the first 
and second terms on the right hand side of Eq. ( |7^ ) 
are negligibly small as evidenced by the large best 
fit values h = 1.56, ^2 = 4.55, and 62 = 0.921. As 
it can be seen in Fig. ^, the addition of the satura- 
tion background (dashed line) only slightly improves 
the systematic difference between the theory and the 
experimental data for \t\ < 4 x 10^^ although it rep- 
resents the data quite well for |t| > 4 x 10~^. 



D. 



XT, Cv, and IA/Ol^^I joint fit 



The good individual fit of the MSR (j)^ model to 
isothermal susceptibility xt, specific heat Cv, and 
coexistence curve \/S.pLy \ has been demonstrated. 
A joint fit of all the three thermal properties leads 
to a complete test of the MSR (fA model with a min- 
imum set of the parameters, {/i, a, Tc, Cb}- Here no 
order parameter saturation was included since its 
correction over the fitting range was small. 

To make sure that no particular measurement 
dominates the joint fit, a proper weighting is needed 
to balance uneven numbers of the experimental data. 
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FIG. 4: Application of the MSR model to the data 
of the ^He coexistence curve. The solid line is the best 
fit with {/^, a} adjusted. The dashed line includes the 
empirical background contribution with {/i, a} fixed from 
the fit without the background. The dot-dashed straight 
line represents the asymptotic prediction from the fit. 



We chose the following weighting in order to normal- 
ize the by the number of data points, 




where N — 



Nc- 



Pl,v 



(75) 



In the joint fit, 



and Cy were fit against temperature T while 
was fit against reduced temperature \t\, and 
/I, a, Cb, and Tc were adjusted. The joint fit results 
are shown in Fig. ^ and Table ||. We note that in 
the joint fit the uncertainties in /i and a are much 
smaller than in the individual fits even though the 
overall goodness of fit is worse in the joint fit. These 
improved uncertainties in ^ and a also lead to the 
improved uncertainties in the critical amplitudes and 
the fluctuation- induced background for specific heat. 

Shown as dashed lines in Fig. || are the Wegner ex- 
pansions to first order with the critical amplitudes, 
F J and Vf , and Af , Bq and Bi , calculated from 
the MSR 0"* model. Bagnuls and Bervillier ||^ have 
argued that the validity range of any cf)'^ model is 
upper-bounded when the difference between the cal- 
culations of the model and the Wegner expansion 
to first order becomes significant. Based on this ar- 
gument, the validity range of the MSR (j)'^ model 
is |t| ~ 1 X 10^^. However, it is interesting to see 
that the MSR (p* model provides a good fit beyond 
|i| = 1 X 10~^ to the experimental measurements of 
the isothermal susceptibility both above and below 
Tc and the specific heat above Tc- 

Close to Tc, the susceptibility data for T > Tc 
and the specific heat data for both T < Tc and 
T > Tc deviate slightly from the theoretical pre- 
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FIG. 5: A joint fit (solid lines) to susceptibility, specific 
heat, and coexistence curve. The fit used all the shown 
Xt data and the data of Cy and |Ap_L,v| over the indi- 
cated range. The dashed lines are the Wegner expansion 
to first order with the listed amplitudes in Table |ll The 
dot-dashed straight lines represent the asymptotic pre- 
dictions from the fit. 



diction. These deviations can be attributed to a 
gravity-induced density stratification. Since the spe- 
cific heat was measured as an average of the whole 
cell while the susceptibility was measured locally 
across a density sensor, there was a stronger grav- 
ity effect in the measured Cy than xt- The grav- 
ity effect on xt{T < Tc) is about a factor of five 
smaller than that on xt{T > Tc) because of the dif- 
ference in XT magnitudes. When Tc is used as an 
adjustable parameter, the individual fits of suscepti- 
bility and specific heat tend to skew Tc such that the 
difference between the experimental measurements 
and theoretical prediction is minimized because of 
the shift in reduced temperatiue for the measure- 
ments. The Tc determined from the fits of the spe- 
cific heat (Fig. m and susceptibility data (Fig. H^) 
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TABLE II: The dimensionless system-dependent param- 
eters for ^He. The adjustable parameters are obtained 
from the joint fit of the 0* model to the measured Xtj 
Cy, and |ApL.v| data of '^He. 



Tc (fit) 
u/u* (fixed) 
X 10* (fit) 

a (fit) 
Cb (fit) 



3.315546 ± 0.000005 
0.999 
1.82 ±0.02 
0.132 ±0.001 
3.74 ±0.02 



^ 



0.146 ±0.001 
0.0296 ± 0.0002 
1.13 ±0.01 
4.95 ± 0.05 



At 
AT 

Bcr 



3.77 ±0.05 
7.05 ± 0.09 
1.01 ±0.01 
0.94 ±0.01 

-5.39 ±0.08 



Bo 



1.019 ±0.002 
0.93 ±0.01 



0.368 ± 0.002 
0.732 ± 0.007 
0.259 ± 0.001 
0.665 ± 0.006 



tends to be higher and lower, respectively, than it 
should be. This tendency was approximately can- 
celled out in the joint fit shown in Fig. |^. The slight 
gravity effect on the experimental measurements for 
1 X 10""* < |t| < 6 X lO"** can be clearly seen in 
Fig.|. 

We mention that in ref. [ p^ earlier measurements 
of the susceptibility of '^He, both above and below Tc, 
were compared with the present data. Also in Table 
I of that reference, the amplitudes of susceptibility 
and coexistence curve data and their ratios, such as 
r^/i?i and r^/r^, obtained from individual fits, 
were presented. 



E. Predictions for correlation length and light 
scattering intensity 

By using u/u*, ^, and a given in Table |l[ the 
dimensionless correlation length can be calculated 
for any given \t\ using Eq. (|5^). Figure ^ shows the 
dimensionless correlation length versus t calculated 
from the MSR (p'^ model for '^He. The length scale to 
recover the dimensional is ^Oi given by Eq. (73). 
Thus one has a dimensional = ^o^o = 2.71 A. This 
value can be directly compared with = 2.6 A mea- 
sured in an acoustic experiment ref. jl^ . Consider- 
ing that the experimental had 10% uncertainty, 
the agreement is very good. 

The correlation length can also be determined 
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FIG. 6: The dimensionless correlation length versus re- 
duced temperature calculated from the MSR 4>^ model 
for ^He. The dot-dashed straight lines represent the 
asymptotic predictions from the fit. 



from a light scattering experiment. Miura, Meyer, 
and Ikushima measured the intensity of scattered 
light of '^He fluid near its critical point ||l^. The in- 
tensity scattered per unit beam length per unit solid 
angle in the fluid, /, is given by 



I = IoAxt sin'^ (jjgik^) , 



(76) 



where Iq is the beam intensity in the scattering 
region, <j> is the angle between the electric field 
of the incident light and the wave vector of the 
scattered light, xt is the susceptibility, and A = 
Tr^kBT{dn^ /dp)'^/XQ. Here n is the index of refrac- 
tion of the fluid, and Aq is the vacuum wavelength of 
the incident light. The function g{k£,) is, for < 10, 
very accurately given by the Ornstein-Zermike ap- 
proximation (1 ± A;^^^)~^+''/^, where k is the scat- 
tering wave vector, ^ is the correlation length, and rj 
is the critical exponent of the fluctuation correlation 
at the critical point. In ref. k — 5.64x 10"* cm~^. 
At t = 1 X 10~^, the value of correlation length can 

1.63x10-* cm, 



-0.63 



^) = 9.2 < 10 was 



be estimated from ^ = 2.7lAi 
hence the condition k£,{t = 1 x 10" 
satisfied for t > 1 x 10~^. Since B = io^sm' cp is es- 
sentially a constant for the experimental condition, 
one can use the knowledge of xt and ^, based on the 
MSR 0'* model, to fit experimental data of the scat- 
tered intensity, with B as an adjustable parameter. 
As it can be seen in Fig. ^ the agreement between 
the experimental data and theoretical calculation is 
reasonably good. 
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FIG. 7: The intensity of light scattered by ^He versus re- 
duced temperature. The theoretical calculated /, using 
the MSR 0* model, is adjusted with a constant ampli- 
tude for loA while u/v* , fi, and a were fixed from the 
values given in Table |l[ 



IV. DISCUSSION 



upon RG theory ||2^. Within this model, the in- 
ternal constants were adjusted such that the critical 
leading amplitude ratios agreed with the values in 
the Fisher and Zinn column in Table |. This CPM 
model was developed to fit simple fluids as well as 
complex fluid systems that exhibit non-monotonic 
crossover behavior. This non-monotonic crossover 
behavior could be described by the CPM approach 
using a finite cut-off wavelength as an additional fit- 
ting parameter. However, in simple fluid systems, 
like ■^He, crossover behavior of different physical 
quantities can be well described within the frame- 
work of the field theoretical (/)^ model without the 
finite cut-off wavenumber. 

NASA supported microgravity flight experiments 
[ pi] , (ittp://miste. jpl.nasa.gov), which are 
under preparation, will take experimental data of 
the susceptibility, specific heat, and coexistence 
curve in the asymptotic region. Combining these 
microgravity measurements in the asymptotic region 
with ground-based measurements in the crossover 
region should permit a rigorous test of the predic- 
tions of recent renormalization theories. 



In this paper we have used parametric expressions 
to calculate the isothermal susceptibility, specific 
heat, coexistence curve, and correlation length along 
the critical isochore and coexistence curve from the 
asymptotic region to the crossover region. All the 
critical leading amplitude ratios were contained in 
the model as listed in Table |. Using only two ad- 
justable parameters in these theoretical expressions 
for the critical contributions, we fit the theory to re- 
cently obtained experimental data for the isothermal 
susceptibility, specific heat, and early experimental 
data of the coexistence curve and light scattering 
intensity. The agreement between the theory and 
experimental measurements is good. 

Further improvements to the minimal renormal- 
ization scheme are desired, especially the five-loop 
Borel resummations throughout the whole range of 
< u < u* . More accurate Borel resummations 
at the fixed point should also lead to improved cal- 
culations of Cc/iiu*) and Cr(u*) so that the resultant 
critical exponents can be compared with other pub- 
lished values (see Appendix ^ . Theoretical insights 
on non-critical contributions are also needed in order 
to formulate more accurate analytical expressions for 
the background contributions. 

While the present minimal subtraction renormal- 
ization model describes quite well to experimental 
measurements along the critical isochore and coex- 
istence curve, it is not as yet a model for a com- 
plete equation-of-state. Recently, Agayan et al. have 
developed a phenomenological crossover parametric 
model (CPM) equation-of-state that is also based 
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APPENDIX A: THE MSR MODEL 
CONSTANTS 

The field-theoretic functions, Cr{u), Ci^(u), and 
Puiu), and the amplitude functions, P+(u), Q{u), 
/±(^^)j F±{u), ^(w), and B{u) are known up to five- 
loop order from expansions around u = 0. However, 
these expansions do not converge. To overcome this 
difficulty, these quantities were expanded to two- 
loop order and then extrapolation terms were added 
to have the functions agree with the calculations of 
high-order Borel resummations at the fixed point j^] . 
All these functions have at least one extrapolation 
term to match the function's value at the fixed point 
u*; some functions also have a second extrapolation 
term in order to match the value of its derivative 
at the fixed point. Listed in this appendix are the 
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values of these extrapolation coefficients, their ori- 
gins, and recent improvements. The effects of these 
coefficient values on the critical exponents and the 
fitting quality in this work are discussed. 

The extrapolation coefficients for the field- 
theoretic functions, Cr{u), C^C'^)? ^-nd (3u{u) in 
Eqs. (|), (0), (|) are ai = 3075, 03 = 30390, 
03 — 37.5, a4 = 14.10, and — 31.85. They are 
taken from Table 2 of ref. . 

The fixed point value for u* ~ 0.040485 is solved 
from the condition = using the given values 

for a4 and 05. The latest published u* value for 
n = 1 is §1 

w* = 0.0404 ±0.0003. (Al) 

The asymptotic critical exponents are linked to the 
exponent functions C^, C,(j,, and /3„ by 



77 = -C0(^^*)--C; = 0.0367, 



(A2) 



v=[2- Cr{u*)]-^ = (2 - C)"' = 0.629 , (A3) 



dPu{u,e = 1) 



du 



= 0.797. 



(A4) 



Once the critical exponents, 77, v, and u are known, 
the remaining important critical exponents can be 
obtained from scaling using 



1 - 2C 

2-C 



2(2-C) 



2 + G 



2-C 



0.112. 



0.326, 



1.235. 



(A5) 



(A6) 



(A7) 



7 = 1.2396 ±0.0013, 



uj = 0.799 ±0.011 



0.504 ±0.008, 



(A13) 



(A14) 



(A15) 



A clear difference exits for the value of the critical 
exponent 7 which warrants further efforts from the 
theoretical community for improvements in the MSR 
model calculation. 

The amplitude function Q(u(l)) for reduced tem- 
perature is expressed as 



Q{u{l))^2 du'-±^exp du- 

Pu[u') Ju{l) 



I3u{u") 
(A16) 

At the fixed point u*, th ere is an identity v ^ — 
2 — C * that simplifies Eq. (A16) and leads to 



Krause et a. 
[Eq. (K.A28)] 



dQ 

du 



Q* 



dP I 



2vPl . (A17) 
obtained the expression 



du 



(AI8) 



with oj — dj3u/du\^, . They also provided a one-loop 
expression for P{u) using a higher-order approxima- 
tion [Eq. (K3.2)] |] 

P+{u) = l-6u{l + bpu). (A19) 

The latest calculation by Larin et al. [Q for n = 1 
gives 

P; = 0.7568 ± 0.0044 . (A20) 

If the theoretically calculated critical exponent i' for 
n — 1, given by Guida and Zinn- Justin [|j, is used 
[H, one has 



A = t/w = 0.502. 



(A8) 



For n = 1, the latest theoretically calculated critical 
exponents given by Guida and Zinn- Justin H are. 



V = 0.6304 ±0.0013, 



T] = 0.0335 ± 0.0025 , 



a = 0.109 ±0.004. 



= 0.3258 ±0.0014, 



(A9) 



(AlO) 



(All) 



(A12) 



b+ = 2j/P_; = 0.9542 ± 0.0059, (A21) 
&_ = 3/2 -2j/P_; = 0.5458 ±0.0059. (A22) 

For the extrapolation coefficients in the expression 
for Q{u), Eq. (p^), Krause et al. ||] deterrnined 
bq = 28.2 and cq = 7.66 such that Eq. ( |aT7| ) 
and ( |AT^ ) were satisfied with the then calculated 
P_^. The values of the extrapolation coefficients 
bq — 20.32 and cq — 6.24 have been readjusted 
to agree with the new Q{u*) — 0.9542. There is no 
value for dP/du\^,, so the new dQ/du\^, has been 
fixed to its old value 

For the amplitude function f±{u) in the expres- 
sion of the susceptibility, Eq. (Eq), b^ = 9.68 comes 
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TABLE III: The values of the various extrapolation coef- 
ficients for the amplitude functions in the MSR model. 



coefficient 


value 


appeared in 


ai 


3075 


Cr{u) for Zr{u) 


a2 


30390 


Cr{u) for Zr{u) 


as 


37.5 


C,piu) for Z,f,{u) 




14.10 


B for 7 7,^(ii\ 


as 


31.85 


(3u{u) for Zu{u)Z4,{u) 


bQ 


20.32 


Q{u) for t{l) 


CQ 


6.24 


Q{u) for til) 


bx 


9.68 


f+(u) for Xt 


dx 


-11.18 


f~ (u) for Xt 


bp 


-5.0726 


F+{u) for C+ 


dp 


-4.6736 


F- lu) for 


be 


-20.6817 


B{u) for Cy 




0.702 


f4,{u) for ApL,v 



from Table 1 of ref.||] and = —11.18 comes from 
Table4ofref. §. 

For the amplitude function in the expression of 



the specific heat, Eqs. (|2D, (l36|), and (gT]), the five- 
loop approximation with a Borel resummation gives 



u*F^ (u*) = 0.3687 ± 0.0040 . 



(A23) 



By combining Eqs. <^)h, (|Al|), and ( |A23D , the old 
interpolation coefficient, dp = —4.04 (Table 4 of 
ref.0), becomes dp = —4.6736. The latest five-loop 
calculation also gives 

u*[F_{u*) - F+{u*)] = 0.4170 ±0.0036. (A24) 



Using Eqs. ( [A23| ) and ( [A24| ), one has 
- u*F+{u*) = 0.0483 ± 0.0076 . 



(A25) 



By combining Eqs. (|37|)a, (|A1|), and ( |A25| ), the 
old interpolation coefficient, bp = 5.04 (Table 1 of 
ref.Q), becomes bp = —5.07. 
We modify Eq. (H) to be 



B{u) = ^ + 9{1 + bBu)u^ 



(A26) 



with bB = —20.68 in order to satisfy the five-loop 
Borel resumed results [Eq. (L2.34)] g 



B{u*) = 0.5024 ±0.001 . 



(A27) 
derived in 



All the calculations use the value of u 
this paper. 

Table ID lists the values of the various extrapo- 



lation coefficients for the amplitude functions in the 
MSR (/)■* model. Table lists the values of the var- 
ious amplitude functions at the fixed point u*. 

Equations (|2^) and ( |26|) provide a clear identifica- 
tion of the leading critical divergence and crossover 



TABLE IV: The values of the various amplitude func- 
tions at the fixed point u* . These values are used in the 
calculation of the leading critical amplitude ratios with 
u* = 0.040485 . 



coefficient 


value 


appeared in 


bl 
bl 


0.9542 ± 0.0059 
0.5458 ± 0.0059 


t{l) in Eq. (20) 
t{l) in Eq. (~) 


/; 
/- 


0.9767 
2.413 


rj for Xt in Eq. S 
Fq for Xt in Eq. (p2|) 


u*Fl 
B* 


0.0483 ± 0.0076 
0.3687 ± 0.0040 
0.5024 ±0.001 


for C; in Eq. ( 
Aq for in Eq. ( 
A* for in Eq. ( 


45|) 

la) 

40^ 




3.175 


Bo for ApL,v in Eq. (|54|) 



contribution in a multiplicative form. In the origi- 
nal expressions, the critical divergence is contained 
implicitly in the integrals of C,r and in Eqs. ( [Tsl ) 
and (|2^). The calculated Cr(u*) and Crpi"^*) using 
Borel resummations at the fixed point lead to the 
critical exponents v and ry that are slightly differ- 
ent from the latest values given by Guida and Zinn- 
Justin. Because of the expressions in Eqs. (|9|) and 
(p6|), the critical exponent values given by Guida and 
Zinn- Justin are used for the leading divergence. The 
inconsistency is only in the crossover part in the in- 
tegrands of [Criu) — Cr{u*)] and [C^{u) — (^^{u*)] that 
go to zero as the fixed point is approached. 



APPENDIX B: DERIVATION OF 
SUSCEPTIBILITY AMPLITUDES 

Expressions for the Wegner expansion of the sus- 
ceptibility will be derived in this Appendix that were 
not presented in previously published work. Multi- 
plying Eq. (|26| ) by Eq. (|20| ) to the power 7 yields 



X±\tr = Xo [b±{l)to] 



exp[-F^(Z)-7j;(0] 
f±{l) 



(Bl) 

In order to expand the exponent functions, Fr{u{l)) 
and F^{u{l)), based on Eqs. (|2^) and (|2^), one needs 
to expand first the function for the flow equation, 
l3u{u{l)), to the first order in [u{l) — u*], 



f3u{u{l)) = cu[u{l) - u*] + 0[{u{l) - u*f] , 

where lo = d(3u/du\u* and /3„(u*) = 0. 
Fr{u*) — and F^(u*) = 0, one obtains 



(B2) 
Since 



F,Ml)) = hm 



Cr{u{l)) ~ CAU*) 



uj[u{l) 



.(0 



(B3) 

uH) ~ u*] + 0[{uil) - u*f] , 
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F^[u[l)) = hm p-— ^ — [u[l)-u 



u{l)->u' Ll![u{1) — U 

o[{u{i)-u*r] 



(B4) 

u{l)-u*]^0[{u{l)-u* f] , 



/±K0) = f±iun + .r±{unHi)-u*] 

+ 0[{u{l) - u* f] , (B5) 



&±(u(/)) = b±{u*) + b'±{u*)[u{l) - u* 
+ OMl)-ur]- 



(B6) 



The expression b±{u{l))^ exp[—F^{u{l)) — 
jFr^u^l))]/ f±{u{l)) is then expanded in terms 
of [u{l) ~ u*], dropping the higher orders, to give 



7 



f±{u' 



(B7) 



uj b± 



^ f± 



[nil) - u*] 



The solution of the flow equation with /?„ (u) approx- 
imated by Eq. (p^) is 



u{l)-u* 



u — 



(B8) 



By expressing I in terms of |i| and dropping higher 
order terms, one has / = \t\'^ / [h'jJ^oY ^'^d 



(B9) 



1-7 77— + — + — 

' b± UJ f± 



u — u , ,, 

X i^r 



Comparing Eq. (|B9|) to the standard Wegner expan- 
sion to the first term, sec Eq. (jsi]), one obtains the 
critical amphtudes of the susceptibihty expressed 
analyticaUy in Eqs. (js^) and (^3|) with A = vlo. 



APPENDIX C: DERIVATION OF SPECIFIC 
HEAT AMPLITUDES 

A derivation of the critical amplitudes and con- 
stant background of the specific heat in the additive 
renormalization form will be given in this Appendix. 
This derivation is consistent with the one for sus- 
ceptibility given above and is different from the one 
given by Schloms and Dohm B . 



First an expansion expression for the function 
A{u{l)) will be derived that is an approximate so- 
lution of Eq. (p8[). By expanding B{u) and Cr{u) 
around u* and omitting higher order terms beyond 
the linear term, Eq. (Isl) becomes 



dl 



(CI) 



+ [a/v~2C{u*){u-u*)]A{l). 



Then Eq. ( |B8| ) is used to replace {u — u*) with Z", 
yielding 



dAil) 
dl 



where 



H + Yl'^ + {G + ZI'^)A{1). 



H^AB{u*) 

V 

Y = AB'{u){u - u) 



With a variable change of 



(C2) 



(C3) 
(C4) 

(C5) 
(C6) 



(C7) 



Eq. (C2) becomes 

H Y 



dA{v) 
dv 



Z 



= - + -v+l-+v]Aiv). (C8) 



The solution of Eq. (|C8D is 



A(v) ~ exp(w) v' 



(C9) 













X 















II 

z 



G 

UJ 



where Ki is a constant to be deter mine d through 
the initial condition. Expanding Eq. ( |CS| ) in v and 
keeping only the linear terms of 1"^'^ and one 
obtains 



A{1) = A{u*) + K2{u - u*) r + Ks l"/"" , 
where Eq. (Eol) is used for A{u*) and 



1 



vY 



A* 



A 



A-: 



{u — u*) 
-^[2B'{u*)~A*C{n*)] 



(CIO) 



(Cll) 



K'i in Eq. (CIO) will be eliminated through initial 
condition at the reference point I — 1 



^3 = Ai- A* - K2{u-u*), 



(C12) 
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where Ai = A{1 = 1). Since Ai has not been given 
as a fitting parameter, it is calculated from the nu- 
merical solution of Eq. (3^) with A{u = 0) = —2. 
Substituting Eq. (|cTo|)into Eq. (H), one has 



Co 



X [F±{l)-A* -K2{u-u*)l'- 



(C13) 



Replacing P/"^ in Eq. ( |cT3| ) with t/to from Eq. (gj 
leads to 



5 
Co 



:-exp[2F,(0]K3 (C14) 

+ exp[(2-a)F,(0] [biiW \t\-" 
X [F±(/) - A* - K2{u - . 



Fr{l), F±{1) , an d ^±(0 are expanded according to 
Eqs. (p3|), (B5), and ( p^ ) respectively. Higher or- 
der terms than 0[(u(/) — u*)^] or 0(l^'^) are dropped 
in the e xpa nsion, and (u{l) — u*) is replaced us- 
ing Eq. (B8). By using the approximation of l'^ = 
\t\^ / (blto)^ , one finally has 



C 

-± = -K3+ (f; - A*) {bltoy |tr" 



X + 



/■I h' 
+ (2-a)^ + «-± 
uj b± 



X {u — u*) 



\t\' 



(C15) 



By comparing Eq. (C15) to the standard Wegner 
expansion to the first term, see Eq. (|4^) , one obtains 
the analytical expressions for the critical arnplitudes 
of the specific heat given in Eqs. (^ ) and (p6|). The 
critical background specific heat is also identified as 



Bcr 



-C0K3 



(C16) 



APPENDIX D: Xt AND Cy EXPERIMENTAL 
MEASUREMENTS 



We list in this Appendix the dimensionless exper- 
imental measurements of isothermal susceptibility, 
specific heat, and coexistence curve of '^He. The 
ITS90 temperature standard was used in the follow- 
ing tables. 



TABLE V: The dimensionless experimental measure- 
ments of the "^He isothermal susceptibility Xt- = 
3.315545K was obtained from the joint fit of Xt ^'Hd 
Cv to the MSR(?i'' model. The index = 0, 19, 8 cor- 
responds respectively to T > Tc,{T < Tc, liquid}, 
{T < Tc, vapor}. 



T T/Tc - 1 



Xt 



index 



3.3157200 
3.3158500 
3.3158770 
3.3159900 
3.3161700 
3.3161770 
3.3163600 
3.3164000 
3.3166570 
3.3173270 
3.3175600 
3.3179100 
3.3179170 
3.3179400 
3.3180170 
3.3185400 
3.3188970 
3.3202400 
3.3221970 
3.3239900 
3.3259870 
3.3295500 
3.3341869 
3.3404500 
3.3487469 
3.3487600 
3.3596699 
3.3745467 
3.3897402 
3.3999300 
3.3999967 
3.4402800 
3.4999362 
3.4999899 
3.6000160 
3.6000502 
3.8000406 
3.3149265 
3.3144595 
3.3137295 
3.3135883 
3.3125410 
3.3122627 
3.3096911 
3.3050875 
3.2969714 
3.2824717 
3.2630551 
3.2566732 
3.2323114 
3.1836186 
3.1064086 
3.3144595 
3.3137295 



5.288e-05 
9.209e-05 
1.002e-04 
1.343e-04 
1.886e-04 
1.907e-04 
2.459e-04 
2.580e-04 
3.355e-04 
5.376e-04 
6.078e-04 
7.134e-04 
7.155e-04 
7.225e-04 
7.457e-04 
9.034e-04 
l.Olle-03 
1.416e-03 
2.006e-03 
2.547e-03 
3.149e-03 
4.224e-03 
5.623e-03 
7.512e-03 
l.OOle-02 
1.002e-02 
1.331e-02 
1.780e-02 
2.238e-02 
2.545e-02 
2.547e-02 
3.762e-02 
5.561e-02 
5.563e-02 
8.580e-02 
8.581e-02 
1.461e-01 
-1.864e-04 
-3.273e-04 
-5.475e-04 
-5.901e-04 
-9.059e-04 
-9.899e-04 
-1.765e-03 
-3.154e-03 
-5.602e-03 
-9.975e-03 
-1.583e-02 
-1.776e-02 
-2.510e-02 
-3.979e-02 
-6.308e-02 
-3.273e-04 
-5.475e-04 



2.763e+04 
1.524e-f04 
1.278e+04 
9.468e+03 
6.006e+03 
5.600e-f03 
4.477e+03 
4.457e+03 
2.966e+03 
1.670e+03 
1.469e+03 
1.256e+03 
1.180e+03 
1.194e+03 
1.120e-H03 
9.472e-^02 
7.790e+02 
5.268e+02 
3.450e+02 
2.587e+02 
1.970e+02 
1.407e+02 
9.850e+01 
7.046e+01 
4.980e+01 
4.979e+01 
3.535e+01 
2.480e+01 
1.884e-f01 
1.642e+01 
1.620e+01 
1.014e+01 
6.540e+00 
6.540e+00 
3.990e+00 
3.939e+00 
2.126e+00 
1.151e+03 
8.291e+02 
3.271e+02 
3.665e+02 
2.145e+02 
1.868e+02 
9.418e+01 
4.708e+01 
2.344e+01 
1.225e+01 
7.408e+00 
6.555e+00 
4.455e+00 
2.696e+00 
1.717e+00 
6.347e-^02 
3.326e+02 








































19 
19 
19 
19 
19 
19 
19 
19 
19 
19 
19 
19 
19 
19 
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TABLE VI: The dimensionless experimental measure- 
ments of the ^He specific heat Cy. Tc = 3.315545K 
was obtained from the joint fit of Xt ^^'^ C'v to the 
MSR (j)* model. 





T 


T/Tc- 


- 1 




3. 


.014570 


-9. 


.078e- 


■02 


7.881 


3. 


.021686 


-8. 


.863c- 


■02 


7.946 


3. 


.028699 


-8. 


.652c- 


■02 


7.959 


3. 


.035582 


-8. 


.444c- 


■02 


8.038 


3. 


.042383 


-8. 


.239e- 


■02 


8.088 


3. 


.049135 
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.035e- 
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.055774 
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.835e- 
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TABLE VII: The experimental measurements of the 
^He reduced density l^pL,v along the coexistence 
curve. The data was provided by Prof. H. Meyer as it 
was used in ref. |l^. The index = 0, 8 corresponds to 
liquid and vapor respectively. 
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